Phase Diagram of Heavy Fermion Metal CeCoIn 5 
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We present a comprehensive analysis of the low temperature experimental H — T phase diagram 
of CeCoIns. The main universal features of the diagram can be explained within the Fermi- liquid 
theory provided that quasiparticles form so called fermion-condensate state. We show that in this 
case the fluctuations accompanying an ordinary quantum critical point are strongly suppressed 
and cannot destroy the quasiparticles. Analyzing the phase diagram and giving predictions, we 
demonstrate that the electronic system of CeCoIn 5 provides a unique opportunity to study the 
relationship between quasiparticles properties and non-Fermi liquid behavior. 
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Although much theoretical efforts have been spared 
to understand the non-Fermi liquid behavior (NFL) of 
heavy fermion (HF) metals using the concept of quantum 
critical points, the problem is still far from its complete 
understanding since the experimental systems display se- 
rious discrepancies with the theoretical predictions Q. 
Common belief is that a quantum critical point (QCP) 
is the point where a second order phase transition oc- 
curs at temperature T — > 0, and where both thermal 
and quantum fluctuations are present destroying quasi- 
particles and generating a new regime around the point 
of instability between two stable phases 0, Q- Recent 
experimental studies of the CeCoIns HF metal provide 
valuable information about the NFL behavior near pos- 
sible QCP due to its excellent tunability by a pressure 
P and/or a magnetic field H 0, 0, 0|. The experimen- 
tal studies have shown that besides a complicated H — T 
phase diagram, the normal and superconducting proper- 
ties around the QCP exhibit various anomalies. One of 
them is power (in both T and H) variation of the re- 
sistivity and heat transport 0, 0, E| Q, , inherent to 
both NFL and Landau Fermi liquid (LFL) regimes. The 
other one is a continuous magnetic field evolution of a 
superconductive phase transition from the second order 
to the first one U El- Above anomalous power laws 
can be hardly accounted for within scenarios based on 
the QCP occurrence with quantum and thermal fluctu- 
ations. For example, the divergence of the normal-state 
thermal expansion coefficient, a/T is stronger than that 
in the 3D itinerant spin-density-wave (SDW) theory, but 
weaker than that in the 2D SDW picture . This brings 
the question of whether the fluctuations are responsible 
for the observed behavior, and if they are not, what kind 
of physics determines the above anomalies? On the other 
hand, the direct observations of quasiparticle band in 
Celrlns have been reported recently [12j. However, if 
the quasiparticles do exist, why they are not suppressed 
by the fluctuations? 

In this letter we show that these problems can be 
resolved within LFL theory provided that quasiparti- 



cles form the so-called fermion-condensate (FC) state 
jl'illl l| emerging behind the fermion condensation quan- 
tum phase transition (FCQPT) We show that near 
FCQPT the fluctuations are strongly suppressed while 
FC by itself is "protected" from above fluctuations by 
the first order phase transition. We analyze the experi- 
mental H — T phase diagram of CeCoIns and show that 
its main universal features can be well understood within 
the theory based on FCQPT. We demonstrate that the 
electronic system of CeCoIns can be shifted from the or- 
dered to disordered side of FCQPT by a magnetic field, 
therefore giving a unique possibility to study the rela- 
tionship between quasiparticles and NFL behavior. 

To study the low temperature universal features of HF 
metals, we use the notion of HF liquid in order to avoid 
the complications related to the crystalline anisotropy of 
solids. This is possible since we consider the (universal) 
behavior related to the power-law divergences of observ- 
ables like the effective mass, thermal expansion coeffi- 
cient etc. These divergences are determined by small (as 
compared to those from unit cell of a corresponding recip- 
rocal lattice) momenta transfer so that the contribution 
from larger momenta can be safely ignored. 

Let us consider HF liquid characterized by the effective 
mass M*. Upon applying the well-known equation, wc 
can relate M* to the bare electron mass M yjj 



M* 



1-N F 1 { P f,Pf)/S' 



(1) 



Here Nq is the density of states of a free electron gas, 
Pi? is Fermi momentum, and F l {pp,pp) is the p-wave 
component of Landau interaction amplitude. Since LFL 
theory implies the quasiparticle density in the form x = 
p F /3ir 2 , we can rewrite the amplitude as F 1 (j>p,pp) = 
F 1 (x). When at some x — x-pci the F 1 (x) achieves some 
critical value, the denominator in Eq. Q tends to zero 
so that the effective mass diverges at T = 0. Beyond 
the critical point Xfc the denominator becomes negative 
making the effective mass negative. To avoid physically 
meaningless states with M* < 0, the system undergoes 
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FCQPT with FC formation in the critical point x — xpc- 
Therefore, behind the critical point xfc the quasiparticle 
spectrum is flat, e(p) = in some region pi < p < Pf 
of momenta, while the corresponding occupation number 
no(p) varies continuously from 1 to 0, < n (p) < 1 [l3j . 
Here fj, is a chemical potential. 

To investigate the FC state at T — 0, we apply weak 



BCS-like interaction [18| with the coupling constant A 
and see what happens with the superconducting order 
parameter k(p) as A — > 0. To do so, we express n(p) via 
the frequency integral from Gorkov function F + . To find 
latter function, we write the usual pair of equations for 
functions F + and G (see e.g. Ref. [l6|). The solution of 
above equations yields 



F + {k,uj) 
G(k,u) = 



i\ < ft(p) > 
27r£i(p) 

l-n(p) 



1 



1 



J5i(p)+*0 u + Ei(p) 
»(P) t.2 



cj — Ei(p) + iO lu + E 1 (p) - iO ' 



^( P ) = ( e (p)-//r+A^ 



(2) 



where the superconducting gap, 

K(P) (2^' (3) 

with 

/OO J 
^F + (p,^)^, (4) 

is related to the dispersion e(p) as 

£ (P)-M = A( P )^H2. (5) 

Here n(p) = v 2 (p), the superconducting order parameter 

«(p) = «(p)«(p) = \A*(p)(l ~ n (P))> U (P) and W (P) are 
the coefficients of corresponding Bogoliubov transforma- 
tion, u 2 (p) + v 2 (p) — 1. Next we observe from Eqs. (3) 
and (5) that when A — > the dispersion e(p) becomes 
flat, while Eq. (2) in the lowest order in A becomes 



F+(k, U ) = -Vno(p)(l-n (p)) - — ^ , 

G (k,u)) = — — + -, e(p)-n = 0. (6) 

LU + IV UJ — IV 

To check this answer, we integrate Fq over frequencies, 

/OO J 
F + (p, — = Vno(p)(l-n (p)), (7) 

and conclude that «(p) = \/no(p)(l ~ "o( P ))- 

Since in our model the transition temperature T c ~ 
A ~ A — ► 0, the order parameter k(p) vanishes at any 
finite temperature so that the quasiparticle occupation 
number is given by Fermi-Dirac function which we rep- 
resent in the form 

•toTj-^D-Ttal^p. ( 8 ) 



Observing that at T — > the distribution function satis- 
fies the inequality < 7io(p,T) < 1 at < p < p/, we 
conclude that both the entropy S and the logarithm in 
the right hand side of Eq. © are finite even at T — > 0. 
In this case, the entropy Snfl is related to the special so- 
lution n (p) and contains the temperature independent 
term So = S^fl(T — > 0). For this special solution n (p) 
the dispersion is flat, e(p) = \x |20f. Thus, the occu- 
pation number no(p) represents the special solutions of 
both BCS and LFL equations determining the NFL be- 
havior of HF liquid. Namely, as it follows from Eq. JJJ), 
contrary to conventional BSC case, the FC solutions are 
characterized by infinitesimal value of superconducting 
gap, A — > 0, while the order parameter /to(p) remains fi- 
nite and the entropy S = 0. At the same time, in contrast 
to the standard solutions of the LFL theory, the special 
ones are characterized by the finite superconducting or- 
der parameter n(p) at T = 0. At T — > both the normal 
state of the HF liquid with the finite entropy S'nfl and 
the BCS state with S — coexist being separated by 
the first order phase transition where the entropy under- 
goes a finite jump SS = So- Due to the thermodynamic 
inequality, 

SQ < T5S, (9) 

the heat SQ of the transition is equal to zero making the 
other thermodynamic functions continuous. Thus, both 
at the FCQPT point and behind it there are no criti- 
cal fluctuations accompanying second order phase tran- 
sitions and suppressing the quasiparticles. As a result, 
the quasiparticles survive and define the thermodynamic 
properties of the HF system. 

On the basis of above special solutions peculiarities 
we can explain the main universal properties of the 
H — T phase diagram of the HF metal CeCoIns shown 
in Fig. 1. The latter substance is a superconductor with 
T c = 2.3 K, while a field tuned QCP with a critical field of 
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H c q = 5.1 T coincides with H C 2, the upper critical field 
where superconductivity vanishes 0, ■ We note that 
in some cases H c0 = 0. For example, CeRu 2 Si 2 shows 
no magnetic ordering down to lowest temperatures |2l| . 
Therefore, in our simple HF model H c q can be treated as 
fitting parameter. It follows from the above considera- 
tion that H c q ~ H C 2 is an accidental coincidence. Indeed, 
H C 2 is determined by A which in turn is given by the cou- 
pling of electrons with magnetic, phonon, etc excitations 
rather than by H c q. As a result, under application of 
a pressure which influences differently A and H c q, the 
above coincidence will be removed. This is in agreement 
with the facts from Ref. Q. 




FIG. 1: H — T phase diagram of CeCoIns. Superconducting- 
normal phase boundary is shown by the solid and dashed 
lines with the solid square showing the point where the super- 
conducting phase transition changes from the second to the 
first order. The solid squares (solid and dashed lines repre- 
sent the calculated dependence of the T*(H) as described in 
the text) obtained from resistivity measurements denote the 
crossover boundaries between the NFL and LFL regimes 
The inset shows the magnetic field dependence of T 2 LFL co- 
efficients of charge A{H) oc H^ 4/3 and heat B{H) <x H~ 4/3 
transport with experimental data taken from Ref. Q. 

At relatively high temperatures, the superconducting- 
normal phase transition in CeCoIns shown by the solid 
line in Fig. 1 is of the second order 0,0. In that case, 
the entropy and the other thermodynamic quantities are 
continuous at the transition temperature T C (H). Since 
H C 2 a H c q, upon the application of magnetic field, the 
HF metal transits to its NFL state down to lowest tem- 
peratures as seen from Fig. 1. As long as the phase 
transition is of the second order, the entropy of SC phase 
Ssc(T) coincides with the entropy <Snfl(T) of NFL state, 

S SC (T^T C (H)) = S NFL (T^T C {H)). (10) 

Since Ssc(T — > 0) — > 0, Eq. JTJJl cannot be satisfied 
at sufficiently low temperatures due to the presence of 
temperature-independent term Sq. Thus, in accordance 



with experimental results |3, |lfjj , the second order phase 
transition converts to first order one below some temper- 
ature T (H). The prediction that the superconducting 
phase transition may change its order had been made 
in the early 1960-s |l9j. This prediction is corroborated 
by our general analysis based on Eq. (|1L)|) . Namely, if 
the superconducting phase were replaced by some other 
ordered phase separated from the NFL phase by the sec- 
ond order phase transition at H = 0, then at some tem- 
perature Tq(H) this phase transition should change its 
order. In our case the NFL phase plays a role of disor- 
dered one. But, as it follows from above consideration, 
the NFL phase has the temperature independent entropy 
term Sq. Since in the ordered phase the Nernst theorem 
(S — > as T — > 0) should hold, we come to conclusion 
that there is the entropy step (from So to zero) as T — > 
while a system traverses the phase transition line from 
ordered phase to NFL one. This means that this phase 
transition should change its order at Tq(H). For exam- 
ple, the AFM phase transition with Tjv (H) (representing 
the field dependence of Neel temperature) should become 
first order at some finite temperature To (if). Besides 
the jump in the entropy, this first order transition, for 
example, should result in a jump in the sample length, 
corresponding to a divergence of a oc —dS/dP. Under 
constant entropy (adiabatic) conditions, there should be 
a temperature step as a magnetic field crosses the above 
phase boundary due to inequality (jSJ). Indeed, the en- 
tropy jump would release the heat, but since S — const 
the heat is absorbed, causing the temperature to decrease 
in order to keep the constant entropy of the NFL state. 
Note that the minimal jump is given by the temperature- 
independent term So, which can be quite large so that 
the corresponding HF metal can be used as an effective 
cooler at low temperatures. 

The entropy Snfl determines the anomalous behav- 
ior of CeCoIns in the NFL region of the phase diagram. 
The term So ~ (j?f —Pi)/pF can be determined from the 
experimental data on spin susceptibility (following Curie 
law) and the specific heat jump AC at T c In HF 

metals like CeCoIns the normalized jump AC/C n ~ 4.5 
is substantially higher than the ordinary BCS value [25| , 
where C n is the specific heat of a normal state. In 
the case of FC, the specific heat jump is not propor- 
tional to T c and is related to the fermion condensate 
parameter Sp^c = ipf ~ P4.)/pf, therefore the normal- 
ized jump AC/Cn can be large |23,|2J|. This estimation 
gives Spfc — 0.044 [glj. The entropy Snfl determines 
also both thermal expansion coefficient a = —dS/dP 
and Griineisen ratio T = a/C n of Fermi liquids with 

fc mum. Since the entropy has the temperature 

independent part So, the thermal expansion coefficient 
a ~ —dSo/dP becomes temperature independent at low 
temperatures. Therefore, at T — > a(T) — > const while 
the specific heat C„(T) -> 0. As a result, T(T -> 0) 
diverges in coincidence with the facts from Ref. [llj| . 
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Now we consider the LFL behavior tuned by a mag- 
netic field H > H c q. Since the NFL behavior of 
CeCoIn 5 coincides with that of YbRh 2 (Sio.95Ge .o5)2 and 
YbRh 2 Si 2 HE 113 we would expect that the LFL behav- 
ior of these substances would also coincide. For example, 
in YbRh 2 (Sio.95Geo.o5)2 the scattering coefficient A(H) 
in the resistivity p — po + AT 2 , with po being the temper- 
ature independent part, diverges as A{H) oc (H — -ff c o) _1 
while in CeCoIn 5 it diverges as A(H) oc (H - H c0 ) c 
with the exponent c ~ —4/3 0,0]. In magnetic fields, 
the exponent c = — 1 characterizes the function A(H) of 
HF liquid with FC [H, while the exponent c = -4/3 de- 
scribes the function A(H) of HF liquid on the disordered 
side of FCjliEi- 

To understand this striking change in the behavior of 
CeCoIns, we recall that FC has just appear in this sub- 
stance since Spfc = {Pf — Pi)/PF — 0.044 <c 1. As 
soon as magnetic field is sufficiently high, H > H CI , (H CI 
is a critical field destroying FC state), Zeeman splitting 
§Pf = (pfi —PF2)/pf of the two Fermi surfaces of HF liq- 
uid exceeds the condensate parameter, Spf > Pfc, an d 
the HF liquid with FC becomes LFL placed on the dis- 
ordered side near QCP. Here pfi and pf2 are the Fermi 
momenta of the two Fermi surfaces formed by the ap- 
plication of a magnetic field. The splitting can be es- 
timated as p f 5pf/M*(H) ~ Hp,B, where p-B is Bohr 
magneton. Taking into account that A(H) oc (M*(H)) 2 
we obtain (H cr — H c q)/H c o ~ (ciSpp) 3 - Our estimations 
of the coefficient c\ based on the experimental function 
A(H) show that c\ ~ 5, and we obtain that reduced 
field {H cr - H c0 )/H cQ ~ (cxdpp) 3 ~ 0.02. Thus, we 
can safely suggest that the reduced field of 0.02 is much 
smaller then minimal reduced field 0.1 where A(H) mea- 
surements have been carried out in Ref. 0, 0- As a 
result, the electronic system of CeCoIn 5 is placed on the 
disordered side of FCQPT by the application of such high 
field and reveals A(H) oc (H — iJ c o)~ 4 / 3 . We can see 
from the inset to Fig.l that the coefficient B(H) has the 
same critical field dependence. Here B(H) stands for 
the Independent contribution to the thermal resistiv- 
ity and related to A{H) by a field-independent factor, 
A(H) I B(H) ~ 0.47, as it should be in the case of ordi- 
nary metals 0. We conclude that the LFL behavior of 
CeCoIns corresponds to the LFL behavior of HF liquid 
placed on the disordered side of FCQPT |3l| • At decreas- 
ing field when H < H cri we predict that the exponent c 
will change from c = —4/3 to c = —1. 

At finite temperatures, the system remains in the LFL 
state, but there exists a temperature T* (H) where the in- 
fluence of FC is recovered and the related NFL behavior 
is restored. To calculate the function T*(H), we observe 
that the effective mass M* cannot be changed at T*(H). 
Since at T > T*{H) the effective mass M*{T) oc 1/T 
and at T < T*(H), M*{H) oc (H - H c0 )- 2 / 3 , we have 
T*{H) oc (H — H c0 )~ 2 / 3 . In high magnetic fields, there 
is a new crossover line because the effective mass starts 



to depend on temperature as M*(T) cc T -2 / 3 [H H3 
and T*{H) becomes T*(H) oc (H - H c0 ). As it is seen 
from Fig. 1, the obtained results are in good agreement 
with experimental facts. It is worth to note that the re- 
covery of FC state can be observed in measurements of 
a tunneling conductivity which is expected to be notice- 
ably asymmetrical with respect to the change of voltage 
bias from V to — V in Fermi systems with FC This 
asymmetry can be observed when the HF metal with FC 
is either normal or superconducting because the asymme- 
try is indeed determined by the violation of the so called 
particle-hole symmetry, while this symmetry holds within 
the LFL theory. This symmetry is not supported by FC 
occupation number no(p, T) since it does not evolve from 
the Fermi-Dirac distribution function. Therefore, the 
tunneling conductivity is asymmetrical when HF metal 
is superconducting, it does not change at T = T c and 
gradually vanishes at rising temperatures. We note also 
that such behavior of the conductivity has recently been 
observed experimentally in CeCoIns On the other 

hand, the behavior of the conductivity can be different 
when the HF metal transits from its LFL to NFL states. 
We predict that in the case of CeCoIn 5 the conductivity 
being symmetrical in the LFL regime becomes gradually 
asymmetrical reaching its maximum in the NFL state at 
elevated temperatures when T > T*(H) and eventually 
vanishes. 

We thank P. Coleman for stimulating discussions. This 
work was supported in part by RFBR, project No. 05- 
02-16085. The visit of VRS to Clark Atlanta University 
has been upported by NSF through a grant to CTSPS. 



* Electronic address: vrshag@thd.pnpi.spb.ru 
t URL: http://cs.uni.opole.pl/~stef Electronic ad- 
dress: stef@math.uni.opole.pl 
[1] G.R. Stewart, Rev. Mod. Phys. 73, 797 (2001). 
[2] S. Sachdev, Quantum Phase Transitions (Cambridge, 

Cambridge University Press, 1999). 
[3] P. Coleman and A.J. Schofield, Nature 433, 226 (2005). 
[4] J. Paglione et al, Phys. Rev. Lett. 91, 2 46405 (2003); J. 

Paglione et al., unpublished (cond- mat/0405157). 
[5] A. Bianchi et al, Phys. Rev. Lett. 91, 257001 (2003). 
[6] F. Ronning et al., unpublished (cond-mat/0602089 ). 
[7] A. Malinowski et al, Phys. Rev. B 72, 184506 (2005). 
[8] E. D. Bauer et al, Phys. Rev. Lett. 94, 047001 (2005). 
[9] K. Izawa et al, Phys. Rev. Lett. 87, 057002 (2001). 
[10] A. Bianchi et al, Phys. Rev. Lett. 89, 137002 (2002). 
[11] N. Oeschler et al, Phys. Rev. Lett. 91, 076402 (2003). 
[12] S. Fujimori et al, unpublished (cond-mat/0602296). 
[13] V.A. Khodel and V.R. Shaginyan, JETP Lett. 51, 553 
(1990). 

[14] G. E. Volovik, JETP Lett. 53, 222 (1991). 
[15] M.Ya. Amusia and V.R. Shaginyan, Phys. Rev. B 63, 
224507 (2001). 

[16] E.M.Lifshitz and L.P.Pitaevskii, Statistical Physics (Part 
2, Butterworth-Heinemann, Oxford, 1999). 



■5 



[17] M. Pfitzner and P. Wolfle, Phys. Rev. B 33, 2003 (1986). 
[18] J.Bardeen, L.Cooper, J.R.Schrieffer, Phys. Rev. 108, 
1175 (1957). 

[19] K. Maki and T. Tsuneto, Progress of Theretical Physics 

31, 945 (1964); K. Maki, Phys. Rev. 148, 362 (1966). 
[20] M. E. Zverev, JETP Lett. 65, 863 (1997). 
[21] D. Takahashi et al, Phys. Rev. B 67, 180407 (2003). 
[22] M.Ya. Amusia, A.Z. Msezane, and V.R. Shaginyan, Phys. 

Lett. A 320, 459 (2004). 
[23] VA. Khodel, M.V. Zverev, and V.M. Yakovenko, Phys. 

Rev. Lett. 95, 236402 (2005). 
[24] M.Ya. Amusia, S.A. Artamonov, and V.R. Shaginyan, 

JETP Lett. 74, 435 (2001). 



[25] C. Petrovic et al, J. Phys. Condens. Matter 13, L337 
(2001). 

[26] J. Custers et al, Nature 424, 524 (2003). 

[27] P. Gegenwart et al, Phys. Rev. Lett. 94, 076402 (2005). 

[28] V.R. Shaginyan, JETP Lett. 79, 286 (2004). 

[29] J. W. Clark, V. A. Khodel, and M. V. Zverev, Phys. Rev. 

B 71, 012401 (2005). 
[30] P. Nozieres, J. Phys. I (France) 2, 443 (1992). 
[31] V.R. Shaginyan, JETP Lett. 80, 263 (2004). 
[32] V.R. Shaginyan, JETP Lett. 81, 222 (2005). 
[33] W. K. Park et al, Phys. Rev. B 72, 052509 (2005). 



